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Abstract:  This  paper  is  focused  on  analyzing  the  impact  of  a  two- 
dimensional  metal  hole  array  structure  integrated  to  the  back-illuminated 
quantum  dots-in-a-well  (DWELL)  infrared  photodetectors.  The  metal  hole 
array  consisting  of  subwavelength-circular  holes  penetrating  gold  layer 
(2D-Au-CHA)  provides  the  enhanced  responsivity  of  DWELL  infrared 
photodetector  at  certain  wavelengths.  The  performance  of  2D-Au-CHA  is 
investigated  by  calculating  the  absorption  of  active  layer  in  the  DWELL 
structure  using  a  finite  integration  technique.  Simulation  results  show  that 
the  performance  of  the  DWELL  focal  plane  array  (LPA)  is  improved  by 
enhancing  the  coupling  to  active  layer  via  local  field  engineering  resulting 
from  a  surface  plasmon  polariton  mode  and  a  guided  Labry-Perot  mode. 
Simulation  method  accomplished  in  this  paper  provides  a  generalized 
approach  to  optimize  the  design  of  any  type  of  couplers  integrated  to 
infrared  photodetectors.  Experimental  results  demonstrate  the  enhanced 
signal-to-noise  ratio  by  the  2D-Au-CHA  integrated  LPA  as  compared  to  the 
DWELL  LPA.  A  comparison  between  the  experiment  and  the  simulation 
shows  a  good  agreement. 
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1.  Introduction 

In  the  past  decade,  the  infrared  detectors  based  on  intersubband  transitions  in  quantum  dots 
(QDs)  have  attracted  much  attention  due  to  lower  dark  currents  and  increased  lifetimes 
compared  with  quantum  well  infrared  photodetectors,  which  are  in  turn  due  to  three- 
dimensional  confinement  and  the  reduction  of  scattering,  respectively  [1,2].  These  epitaxially 
self-assembled  and  nanometer-scale  QD  islands  form  spontaneously  on  a  semiconductor 
substrate  because  of  lattice  mismatch  [3].  More  recently,  infrared  detectors  based  on  QDs 
embedded  in  a  quantum  well  (QW),  known  as  dots-in-a-well  (DWELL)  heterostructure,  have 
been  proposed  [4].  The  DWELL  architecture  brings  the  advantages  of  QWs,  such  as  control 
over  the  operating  wavelength,  while  maintaining  the  advantages  of  QDs,  such  as  normal 
incidence  operation  [5].  In  parallel,  focal  plane  array  (FPA)  development  for  infrared  imaging 
has  proceeded  from  the  first  to  third  generations  (linear  arrays,  2D  arrays  for  staring  systems, 
and  large  format  with  enhanced  capabilities,  respectively).  For  a  step  closer  to  the  next 
generation  of  FPAs,  we  envision  that  perforated  metal  film  structures  [6,7]  will  improve  the 
performance  of  FPA  by  enhancing  the  coupling  to  photodetectors  via  local  field  engineering, 
and  will  enable  spectral  sensitivity  modifications.  In  regard  to  the  improved  performance  at 
certain  wavelengths,  it  is  worth  pointing  out  the  structural  difference  between  previous 
perforated  metal  film  integrated  front-illuminated  single  pixel  devices  [8-10]  and  our  back- 
illuminated  device  (FPA)  [11].  Apart  from  the  pixel  linear  dimension,  it  is  a  distinct 
difference  that  there  is  a  metal  cladding  (composed  of  a  number  of  metals  for  ohmic  contact 
and  the  read-out  integrated  circuit  hybridization)  in  the  FPA  between  the  heavily  doped 
gallium  arsenide  (GaAs)  used  as  the  contact  layer  and  the  read-out  integrated  circuit  (i.e.,  n- 
type  GaAs/QD-absorber/n-type  GaAs/metal  cladding/ROIC);  on  the  contrary,  the  front- 
illuminated  single  pixel  device  consists  of  two  heavily  doped  contact  layers  separated  by  the 
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QD-absorber  on  a  semi-infinite  GaAs  substrate  (i.e.,  n-type  GaAs/QD-absorber/n-type 
GaAs/substrate). 

In  this  paper,  a  single  metal  layer  is  chosen  as  the  metal  cladding  in  the  simulation  to 
allow  clearer  demonstration  of  the  effects  of  a  two-dimensional  square  array  of  circular  holes 
in  Au  layer  (2D-Au-CHA).  As  a  result  of  the2D-Au-CHA  and  metal  ground  plane,  the  2D- 
Au-CHA  integrated  DWELL  LPA  (2D-Au-CHA: DWELL  LPA)  yields  a  wavelength  selective 
resonant  enhancement  of  the  detector  response  in  the  mid-infrared  spectral  region.  The 
analysis  and  comparison  of  modeling  results  with  experiment  indicates  a  good  quantitative 
agreement  of  the  improved  performance  of  DWELL  FPA  measured  using  the  signal  to  noise 
(SNR)  and  the  noise  equivalent  difference  in  temperature  (NEDT)  at  resonance  wavelengths. 
The  improved  performance  results  from  a  surface  plasmon  polariton  (SPP)  mode  and  a 
guided  Fabry-Perot  mode,  that  enhance  x  or  y  (along  the  polarization  direction  used  in 
simulation)  and  z  (along  the  growth  direction)  components  of  electric  field  in  the  active 
region  compared  with  the  unmodified  DWELL  FPA  at  certain  wavelengths. 

2.  Structure  of  2D-Au-CHA  integrated  DWELL  based  back-illuminated  infrared 
photodetector 

In  order  to  find  the  optimum  dimension  of  2D-Au-CHA  to  match  the  first  order  SPP 
resonance  with  the  maximum  of  DWELL  FPA’s  spectral  response,  we  employed  both  the 
commercial  software  (CST  Microwave  Studio:  Computer  Simulation  Technology  GmbH, 
Darmstadt,  Germany)  [12]  based  on  a  finite  integration  technique  and  a  rigorous  coupled 
wave  analysis  [13].  Comparable  results  were  obtained  for  both  techniques  (not  shown  here). 

Indium  arsenide  (InAs)  QDs  were  grown  with  VG-80  solid  source  molecular  beam 
epitaxy  (MBE)  system  using  a  Stranski-Krastonow  growth  mode  [14].  After  the  MBE 
growth,  the  FPA  was  fabricated  (a  dry  etch  to  form  the  mesa,  surface  passivation  (SiNx), 
ohmic  metal  evaporation,  under  bump  metallization,  indium  deposition,  and  indium  reflow 
process)  and  then  hybridized  to  the  read-out  integrated  circuit  by  flip-chip  bonding 
processing.  Following  the  removal  of  substrate  and  AlAs  etch  stop  layer,  the  2D-Au-CHA 
was  fabricated  on  a  1  pm  thick  contact  layer  of  the  DWELL  FPA  as  shown  in  Fig.  1(a).  The 
schematic  diagram  in  Fig.  1(a)  illustrates  a  pixel  in  the  2D-Au-CHA:DWELL  FPA.  The 
zoomed- in  image  in  Fig.  1(b)  depicts  the  DWELL  hetero structure  within  a  pixel.  The 
DWELL  structure  includes  the  n-doped  GaAs  contact  layers  with  200/1000  nm  thicknesses, 
the  active  region  with  20  stacks  of  InAs  QDs  embedded  in  3.5  nm  thick 
In0.15Ga0.85As/Al0.08Ga0.92As  QWs,  each  separated  by  50  nm  thick  Alo.08Gao.92As  barriers. 
Figure  1(c)  shows  geometrical  parameters  of  a  simplified  unit  cell  of  the  2D-Au- 
CHA:DWELL  FPA  used  in  the  simulation.  The  details  of  MBE  growth/fabrication  of 
DWELL  FPA  and  2D-Au-CHA  are  available  in  Ref  [  1 1 , 1 5] . 
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Fig.  1.  Illustration  of  (a)  a  pixel  in  the  2D-Au-CHA:DWELL  FPA  (after  substrate  and  AlAs 
etch  stop  layer  were  removed)  and  (b)  the  DWELL  heterostructure,  (c)  Schematic  view  of  a 
simplified  unit  cell  of  the  2D-Au-CHA:DWELL  FPA  used  for  numerical  simulation. 
Geometrical  parameters  of  2D-Au-CHA: DWELL  FPA:  p  =  1.8  pm;  d=p/2\  tm  =  50  nm;  tup  =  1 
pm;  tabs  =  1-12  pm;  tdown  =  0.2  pm;  tmgp  =  0.2  pm.  In  addition,  the  configuration  of  polarization 
and  propagation  is  depicted  ( E ,  H,  and  k  denote  electric  field,  magnetic  field,  and  wave  vector, 
respectively). 


3.  Simulated  absorption,  responsivity  and  enhancement  ratio 

For  the  purpose  of  understanding  our  experimental  results  quantitatively  and  qualitatively,  we 
performed  3D  full  field  electromagnetic  simulations  of  2D- Au-CHA: DWELL  FPA  using  a 
finite  integration  technique.  The  boundary  conditions  on  the  faces  (of  a  single  cubic  unit  cell) 
parallel  to  the  propagation  vector  were  set  as  perfect  magnetic/electric  conductor  to 
approximate  a  TEM  plane  wave  propagating  through  2D- Au-CHA :D WELL  FPA  as  shown  in 
Fig.  1(c).  All  geometrical  parameters  are  indicated  in  Fig.  1(c).  For  materials  constituting  2D- 
Au-CHA:  DWELL  FPA,  a  Drude  model  was  applied  for  the  dielectric  parameters  of  Au  [16] 
(1.4  times  larger  scattering  frequency  than  bulk  Au  to  fit  the  experimental  values  obtained 
from  Sandia  National  Laboratory  [17]),  n-type  GaAs  [9,18]  (2  x  1018  cm-3)  and  lossy 
aluminium  [16]  used  as  2D- Au-CHA,  contact  layer  and  metal  ground  plane,  respectively. 
Figures  1(b)  and  1(c)  show  that  20  stacks  of  QD  layers  were  simplified  to  an  effective 
absorber  layer  whose  neff  and  keff  (real  and  imaginary  parts  of  refractive  index  in  effective 
absorber)  were  taken  from  Ref  [9],  measured  quantum  efficiency  [19],  and  spectral 
photoresponse. 

The  absorption  in  the  active  layer  (or  effective  absorber)  is  only  useful  to  determine  how 
much  the  response  (e.g.,  photocurrect,  responsivity,  or  SNR)  is  enhanced  with  2D-Au-CHA, 
i.e.,  the  absorption  of  other  layers  is  parasitic.  The  power  absorbed  per  unit  volume  can  be 
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calculated  from  the  divergence  of  the  Poynting  vector  ( V  •  S )  and  is  equivalent 

to  ^  ,  where  co  is  the  angular  frequency  (2 nc/X),  c  is  the  speed  of  light  in 

vacuum,  s"(X)  is  the  imaginary  part  of  dielectric  permittivity,  and  \E(X)\2  is  the  intensity  of  the 
electric  field  [20].  Therefore,  the  absorption  in  one  of  the  layers  constituting  2D-Au- 

rmXX)\E.af 

CHA:DWELL  FPA  is  calculated  byj — ,v  n  ,v  1  dV.,  where  i  =  2D-Au-CHA,  contact 

layers,  effective  absorber  or  metal  ground  plane.  However,  the  absorption  in  effective 
absorber  should  be  calculated  carefully  because  the  QD  based  infrared  photodetector  shows 
different  response  [3]  depending  on  the  orientation  of  electric  field  intensity  (\EX\2,  \Ey\ 2,  or 
\EZ\ 2).  More  specifically,  the  effective  absorber’s  absorption  due  to  2D-Au-CHA  is 
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where,  \Ejea  is  the  intensity  of  j  component-electric  field  in  the  effective  absorber  and  ~8  is 

used  for  >\  (electric  field  along  the  growth  direction  shown  in  Fig.  1(b)  is  dominant  in 
interacting  with  the  QDs  as  compared  to  the  electric  field  parallel  to  the  perforated  metal 
film)  [21].  Figure  2  shows  the  calculated  absorption  in  each  layer  of  the  considered  2D-Au- 
CHA:DWELL  FPA  (calculated  absorption  of  DWELL  FPA  is  not  shown  here). 
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Fig.  2.  Absorption  calculated  in  the  different  layers  of  the  2D-Au-CHA: DWELL  FPA  as  a 
function  of  the  illumination  wavelength. 

Both  DWELL  FPA  and  2D- Au-CHA: DWELL  FPA  are  illuminated  under  the  blackbody 
at  423  K  (the  temperature  used  for  measurements  [11]).  Using  the  absorption  in  the  effective 
absorber  as  shown  in  Fig.  2,  we  calculate  the  responsivity  ( R ),  Rk(A)  =  AkEA  •  RBB at 423K  ,  where 
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k  is  2D-Au-CHA:DWELL  FPA  (or  DWELL  FPA),  RBB  at  423k  can  be  evaluated  using  Planck’s 
law  for  the  blackbody  at  423  K  and  the  equivalent  spectral  current  density,  J(X)  =  N0(E)  x  e~. 
N0(E)  is  the  incident  photon  flux  and  e~  is  the  electron  charge.  Here,  it  is  assumed  that  an 
electron-hole  pair  is  collected  for  each  photon  absorbed  in  effective  absorber,  which 
corresponds  to  an  internal  quantum  efficiency  equal  to  1.  Figure  3  shows  the  simulated 
responsivities  of  DWELL  FPA  and  2D-Au-CHA:DWELL  FPA.  Moreover,  in  order  to 
estimate  the  influence  of  the  2D-Au-CHA  on  the  enhancement  of  the  DWELL  FPA 
efficiency,  we  illustrate  in  Fig.  3  the  ratio  between  the  responsivity  of  2D-Au-CHA:DWELL 
FPA  over  the  responsivity  of  DWELL  FPA,  which  is  represented  by  the  red  circle  and  agrees 
well  with  the  SNR  measurements  as  will  be  shown  in  Fig.  4.  Clearly  there  are  resonances  at 
-4.56  pm,  -5.45  pm,  and  -6.38  pm.  It  is  of  great  importance  to  investigate  the  electric  field 
distributions  in  the  effective  absorber  at  resonant  wavelengths  because  the  absorption  is 
related  to  the  orientation  and  intensity  of  electric  fields  as  shown  in  Eq.  (1)  and  the  effective 
absorber  is  especially  sensitive  to  Ez  due  to  rj.  In  addition,  understanding  this  mechanism 
helps  us  to  improve  the  design  of  the  perforated  metal  film  integrated  infrared  FPA.  At  4.56 
pm  and  6.38  pm,  Ex  distribution  in  z  direction  is  not  due  to  the  propagating  modes  shown  in 
perforated  metal  film  integrated  front-illuminated  device,  but  rather  to  the  guided  Fabry-Perot 
modes  resulting  from  the  presence  of  metal  ground  plane  in  the  FPA  structure;  Ez  distribution 
in  the  plane  parallel  to  the  perforated  metal  film  exhibits  the  SPP  modes  (1st  and  2nd  order 
SPPs  owing  to  the  interaction  between  the  first  and  second  nearest  apertures,  respectively);  Ez 
in  z  direction  displays  the  exponential  decay  of  a  surface  wave  as  shown  in  the  insets  of  Fig. 
3.  The  responsivity  at  5.45  pm  is  increased  by  the  guided  Fabry-Perot  mode,  evidenced  by  a 
standing  wave  pattern  of  Ex  along  z  direction  depicted  in  the  inset  of  Fig.  3. 


Wavelength  (jim) 


Fig.  3.  Simulated  responsivities  of  the  DWELL  FPA  (gray  square),  2D-Au-CHA: DWELL  FPA 
(black  square)  and  enhancement  ratio  (red  circle)  resulting  from  2D-Au-CHA,  i.e.,  R2d-au- 
cha.dwell  fpa  /  R  dwell  fpa •  Note  that  unity  internal  quantum  efficiency  in  effective  absorbing 
layer  is  used.  The  simulated  electric  field  distributions  at  the  corresponding  resonances 
(resulting  from  surface  plasmon  polariton  and  guided  Fabry-Perot  modes)  are  shown  in  the 
inset  to  the  figure. 
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4.  Signal  to  noise  and  noise  equivalent  difference  in  temperature 

To  accurately  quantify  the  resonant  enhancement,  we  consider  the  SNRs  of  2D-Au- 
CHA:DWELL  FPA  and  DWELL  FPA.  The  signal  is  the  intensity  of  averaged  spectral 
response  and  the  noise  is  the  standard  deviation  of  spectral  response  over  20  frames. 
Subsequently,  the  SNRs  of  2D- Au-CHA: DWELL  FPA  and  DWELL  FPA  were  computed  and 
shown  in  Fig.  4  for  the  entire  wavelength  range.  By  taking  the  ratio  between  the  SNRs  of  2D- 
Au-CHA:DWELL  FPA  and  DWELL  FPA,  we  extracted  three  enhancement  ratios,  1.63  (4.5 
pm),  1.21  (5.5  pm),  and  1.59  (6.1  pm)  by  2D-Au-CHA:DWELL  FPA,  respectively.  The 
overall  agreement  between  experimental  results  and  the  simulated  enhancement  ratios  at  the 
corresponding  resonance  is  apparent  from  Fig.  3.  Error  in  resonance  wavelengths  between 
simulation  and  experiment,  AX  =  (Xsim  -  Xexp)/Xexp,  is  less  than  5%,  which  can  be  attributed  to 
imperfections  and  variation  from  the  fabrication  process.  The  enhancement  ratios  over  the  full 
spectrum  are  shown  in  Fig.  4. 

In  addition,  the  NEDT  for  2D-Au-CHA:DWELL  FPA  at  a  device  temperature  of  77  K 
was  measured  by  changing  the  irradiance  via  a  calibrated  blackbody  source  at  various  scene 
temperatures.  By  definition,  the  NEDT  (in  units  of  mK)  is  a  performance  measure  of  FPA’ s 
sensitivity  which  indicates  the  smallest  uniform  scene  temperature  difference  the  FPA  can 
detect.  In  the  measurement  setup,  we  used  a  narrowband  optical  filter  centered  at  5.95  pm 
with  AX  ~140  nm  in  front  of  the  blackbody  source  for  demonstrating  the  sensitivity 
enhancement  near  the  first  resonance.  NEDT  values  were  obtained  by  using  the  FPA’s  output 
response  ( AV0 )  between  two  temperatures  of  the  blackbody  source  (AT),  along  with  the 
recorded  noise  ( Vn )  at  the  lower  temperature  with  the  following  equation  [22-24]:  NEDT  = 
AT/(AV0/Vn).  We  computed  the  NEDT  averaged  over  121  uniform  and  well-behaved  pixels  on 
the  array.  The  averaged  NEDTs  of  2D-Au-CHA: DWELL  FPA  and  DWELL  FPA  were  77.3 
mK  and  109.1  mK  at  an  irradiance  of  3.65  x  10~5  W/cm2.  By  comparison,  this  corresponds  to 
a  sensitivity  improvement  of  2D-Au-CHA:DWELL  FPA  by  a  factor  of  1.4  times  over 
DWELL  FPA. 
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Fig.  4.  Experimental  SNRs  of  the  DWELL  FPA  (gray  triangle)  and  2D-Au-CHA:DWELL 
FPA  (black  triangle),  and  the  corresponding  enhancement  ratio  (red  diamond)  by  2D-Au- 
CHA:DWELL  FPA  over  the  DWELL  FPA. 

5.  Conclusion 

In  summary,  we  have  investigated  the  effect  of  2D-Au-CHA  on  the  performance  of  an 
infrared  DWELL  FPA  both  numerically  and  experimentally.  The  wavelength  selective 
enhancement  due  to  the  2D-Au-CHA  (or  perforated  metal  film)  in  an  infrared  FPA  could 
have  a  dramatic  impact  on  multispectral  (or  hyperspectral)  imaging  system  construction. 
Understanding  the  underlying  mechanism  by  means  of  a  simulation  may  offer  the  possibility 
to  explore  a  large  variety  of  phenomena  in  the  next  generation  of  infrared  detectors  associated 
with  surface  plasmon  polariton  and  guided  Fabry-Perot  modes,  as  well  as  applications  in 
polarimetric  infrared  imaging  and  the  improvement  of  the  design  of  infrared  detectors 
integrated  with  any  type  of  couplers.  A  useful  result  is  that  the  enhanced  electric  fields  inside 
the  effective  absorber  are  strongly  associated  with  the  thickness  of  contact  layer  and  QD  layer 
(i.e.,  spacer  between  perforated  metal  film  and  metal  ground  plane),  therefore  additional 
design  and  fabrication  refinement  are  needed  for  better  performance  of  the  perforated  metal 
film  integrated  FPA. 
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